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Abstract. We present a particular approach to the non-equilibrium dynamics of 
quantum field theory. This approach is based on the Jaynes-Gibbs principle of the 
maximal entropy and its implementation, throughout the initial-value data, into the 
dynamical equations for Green's functions. We use the 4> 4 theory in the large N limit to 
show how our method works by calculating the pressure for a system which is invariant 
under both spatial and temporal translations. 



INTRODUCTION: JAYNES-GIBBS PRINCIPLE 

The objective of this talk is to present a novel approach to a non-equilibrium dy- 
namics of quantum fields [1] . This approach is based on the Jaynes-Gibbs maximum 
entropy principle [2] , which, in contrast to other approaches in use [3-6] , constructs 
a density matrix p directly from the experimental/theoretical initial-time data (e.g. 
pressure, density of energy, magnetization, ionization rate, etc.). We illustrate our 
method on the 4 theory with the O(N) internal symmetry in the large N limit, 
provided that the non-equilibrium medium in question is translationally invariant. 

To start, we consider the following definition of expectation value of some dy- 
namical operator A: (A) = Tr(pA), where the trace is taken with respect to an 
orthonormal basis of physical states describing the ensemble at some initial time 
t; L . Let us consider the information (or Shannon) entropy S[p] = — Tr(plnp) [2]. 
According to the Jaynes-Gibbs principle, we have to maximize S[p] subject to 
constrains imposed by the expectation value of certain experimental/theoretical 
observables: (-Pj[$, d$\(xi, . . .)) = gi(xi, . . .), i = 1, . . . n , where the operators 
Pi[$,<9$], in contrast to thermal equilibrium, need not to be constants of the mo- 
tion; space-time dependences are allowed. The maximum of the entropy determines 
the density matrix with the least informative content. 
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p=^ ) exp^-±J]ld%X l (x 1 ,...)P^,d^ , (1) 

where \ are the Lagrange multipliers to be determined. The 'partition function' Z 
is Z(\j) = Tr jexp (— Yh=i I Tlj d A Xj\(x\, . . .)Pj[<&, 9<&]) j . In the previous equa- 
tions the time integration is not present at all (i.e. g,i{. . .) are especified at the 
initial time ti). In case when the constraint functions gi{. . .) are known over some 
gathering interval (— r, U) the correspondent integration j\dt should be present 
in p. The Lagrange multipliers \ might be eliminated if one solves n simul- 
taneous equations: gi = —8\a.Z/8\. The solution can be formally written as 

\ = SS G \gi, ...,g n ] \rnax /$9i- 



OFF-EQUILIBRIUM DYNAMICAL EQUATIONS 

In this section we briefly introduce the off-equilibrium dynamical (or Dyson- 
Schwinger) equations. For simplicity we illustrate this on a single scalar field $ 
coupled to an external source J described by the action S"[$] = S[<&\ + / J$. 
Associated with this action we have the functional equation of motion [1,5]: 
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Z[J] = -.J a , (2) 



with Z[J\ = Tr{pT c exp(i J c d 4 xJ(x)&(x)} being the generating functional of 
Green's functions. Here C is the Keldysh-Schwinger contour which runs along 
the real axis from ti to tf (tf > ti, tf is arbitrary) and then back to t«. In (2) we 
have associated with the upper branch of C the index " + " and with the lower one 
the index " — " (in the text we shall denote the indices +/— by Greek letters a, (3). 

Let us define the classical field <p a as the expectation value of the field operator 
in the presence of J: i.e. <p a = (3> Q ). Defining the generating functional of the 
connected Green's functions as Z[J] = exp(iW[J]), the two-point Green's function 
is G a p(x,y) = - 5Ja ^ My) = -i{T c <S>{x)<S>{y)) +*($(*)>(<%)>. Eq.(2) is the first 
one of an infinite hierarchy of equations for Green functions. Further equations 
can be obtained from (2) by taking successive variations with respect to J. True 
dynamical equations are then obtained if one substitutes the physical condition 
J = into equations obtained. 

To reflect the effects of the density matrix in the Dyson-Schwinger equations 
it is necessary to construct the corresponding boundary conditions. 3 Using the 
cyclic property of the trace together with the Baker- Campbell-Hausdorff relation: 
e A Be- A = En=o ^Cn, (where C = B and C n = [A, C n _i]), and setting A = ln(p) 
and B = 3>(xi) with x w = ti we obtain the generalized KMS conditions: 



3 ) Let us remind that at equilibrium the corresponding boundary conditions are the Kubo-Martin- 
Schwinger (KMS) conditions. 
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namely for the two-point Green function we have G + -(x,y) = G- + (x,y) + 
Y^k~i ~t\Ti{pQ(x)Ck(x)} . As an example of the latter relation we can choose the 
particular situation when p = exp(—/3H)/Z, in which case we get the well known 
KMS condition: G+_(x; t, y; 0) = G_+(x; t - i/3, y; 0). 



EXAMPLE: OUT-OF-EQUILIBRIUM PRESSURE 

In order to apply our previous results let us consider the 4 theory with the O(N) 
internal symmetry in the large N limit (also the Hartree-Fock approximation). It 
is well known that, in this limit only two-point Green's functions are relevant 
[1,6,7]. The Dyson-Schwinger equations for G a/ 3 are automatically truncated and 
reduce to the Kadanoff-Baym equations [4]: (□ +h1q + ^G aa (x,x)^J G a/3 (x,y) = 
— 5(x — y){(Jz) a p , where a% is the Pauli matrix; Ao and m are, respectively, the bare 
coupling and the bare mass of the theory. If the system is translationally invariant 
the Fourier transform solves the Kadanoff-Baym equations and the correspond- 
ingfundamental solution reads: G a p(k) = fc 2 + ^2+\"( CT3 ) — 2iri5(k 2 + M. 2 )f a f3(k) , 

where the (finite) M is M 2 = ml+i^-G ++ (0). Function f a p(k) must be determined 
through the generalized KMS conditions. Let us now choose the constraint to be 
used. Keeping in mind that we are interested in a system which is invariant under 
both spatial and temporal translations, we choose the constraint <?(k) = (7Y(k)), 
where 7i = uj^a) (k)a(k), with u>k = Vk 2 + M 2 (notice that in the large iV limit the 
Hamiltonian is always quadratic in the fields). The corresponding density matrix 
then reads 

' -0(k)W(k) , (3) 




(2tt) 3 2^," v ' v ') 

with ( 2 ^32cj fc being the Lagrange multiplier to be determined. According to the 
maximum entropy principle we find that /5(k) fulfils equation 

g(k)= (27T)3^K-l ' (4) 

where V denotes the volume of the system. Eq.(4) can be interpreted as the density 
of energy per mode. Similarly as in the case of thermal equilibrium, /3(k) could 
be interpreted as "temperature" with the proviso that different modes have now 
different "temperatures" . 

The generalised KMS conditions in this case are G + _(k) = e~ l3 ^ h °G_ + (k), and 
so the corresponding / ++ reads: / ++ = [exp(/?(k)c<j fc ) — Let us now con- 

sider a particular system in which g(k) = -^p- exp(c<jfc/(T). In this case a is the 
physical parameter which, as we shall see below, can be interpreted as a "tem- 
perature" parameter. This particular choice corresponds to a system where the 



lowest frequency modes depart from equilibrium, while the high energy ones obey 
the Bose-Einstein distribution (typical situation in many non-equilibrium media, 
e.g. plasma heated up by ultrasound waves, hot fusion or ionosphere ionised by 
sun). In terms of the parameter a the Lagrange multiplier may be written as 

/3(k) = \ + -^J2n=i ( ~T +1 exp(-ncc> fc /o-). Notice that when cu k > a, (3 ~ cr _1 , 
and we may see that / ++ approaches to the Bose-Einstein distribution with tem- 
perature a. However, when Uk ~ o the latter interpretation fails. Instead of the 
parameter a, it may be useful to work with the expectation value of /3(k): 



= fd*k[3(k)e-^ = 1 
{P) fd 3 ke-"k/° a 



MK 2 (M/cr) 



(5) 



where K n is the Bessel function of imaginary argument of order n. An interesting 
feature of Eq.(5) is that it is actually insensitive to the value of M. which is impor- 
tant if one wants to use l/{(3) as a "temperature". The actual behaviour of ((3) is 
depicted in Fig.l 
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FIGURE 1. In (a) we plot Eq.(5) at a — lOOMev, while in (b) we plot the difference of equilib- 
rium and non- equilibrium pressures for m r = lOOMeV. 



Let us now consider the renormalized expression for the expectation value of the 
energy momentum tensor [1]: 

(V)ren = nJ ^M^-h+W - G(k)\ - J ^- d [G ++ (k) + G(k)\ , 

with G being the usual (T = 0) causal Green function and Sm 2 = M 2 — m 2 with 
m r being the (T = 0) renormalized mass. The pressure per particle, in the high 
"temperature" expansion (i.e. for large a or small (f3)) for the system described 
by the density matrix (3) may be worked out either in terms of a, using the Mellin 
transform technique [1,8]: 

r>, \ 1 mix ^ ° 2M2 K(v 2 M 2 <r 4 3\ ~ /, , . , , v , 2 \ 



or in terms of 1/ (f3) using the Pade approximation [1]: 

P((P)) = 0.0681122 (/3)" 4 - 0.0415368 (f3y 2 M 2 + A r (-0.000647 ({3}~ 4 
+ 0.0000164 ((3)~ 2 M 2 ) + O (M 2 \n(M([3)); \ 2 r ) . 

It is interesting to compare the previous two results with the high-temperature 
expansion of the same system in thermal equilibrium [7]: 

= T^_T^ TM^ A, (T*_ _ T^M T 2 M 2 \ / /A1_\\ 
1 J 90 24 12tt 8 \UA 24vr 16tt 2 J V \TAirJJ' 

Particularly, the leading "temperature" coefficients in the first two expansions ap- 
proximate to a very good accuracy the usual Stefan-Boltzmann constant for scalar 
theory. The latter vindicates the interpretation of a and 1/ (f3) as temperatures for 
high energy modes. The behaviour of both P(T) and P(cr) are shown in Fig.l. 

SUMMARY AND OUTLOOK 

One of the main advantages of the Jaynes-Gibbs construction is that one starts 
with constraints imposed by experiment/theory. The constraints directly deter- 
mine the density matrix with the least informative content (the least prejudiced 
density matrix which is compatible with all information one has about the sys- 
tem) and consequently the generalized KMS conditions for the Dyson-Schwinger 
equations. We applied our method on a toy model system (O(N) A0 4 theory), 
in the translationally invariant medium. The method presented, however, has a 
natural potential to be extensible to more general systems. Particularly to media 
where the translational invariance is lost. As an example we can mention systems 
which are in local thermal equilibrium. For such systems it is well known [2,8] that 
equilibrium (5 must be replaced by /3(x) (i.e. temperature which slowly varies with 
position). Obviously one may receive this result from the outlined Jaynes-Gibbs 
principle almost for free. Work on more complex systems is now in progress. 
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